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We present a microscopic theory of fluctuation-mediated pairing mechanism in a two-dimensional 
extended Hubbard model. In contrast to conventional wisdom, odd-frequency spin-triplet pair- 
ing can be stabilized near the spin-density-wave critical point. Favorable conditions for the odd- 
frequency pairing are the presence of geometrical frustration and off-site Coulomb interaction. 
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Exploring unconventional pairing functions of super- 
conductivity in strongly correlated electron systems has 
been an important issue in condensed matter physics. 
The pairing function of two electrons must change sign 
under the exchange of two electrons in accordance with 
Fermi-Dirac statistics. Usually the pairing function is 
an even (odd) function of the momentum when the spin 
state of the pair is a spin-singlet (-triplet), since we pos- 
tulate that the pairing function is even as a function of 
frequency. However, the pairing function can be odd, 
as a function of the frequency where the resulting pair- 
ing function has odd (even) parity in momentum space 
for the spin-singlet (-triplet) pairing. For example, the 
"spin-triplet s-wave" or the "spin-singlet p-wave" with 
"odd frequency" can emerge. 

Odd frequency superconductivity has been discussed 
by Berezinskii for the s-wave spin-triplet case[2J in the 
context of 3 He. Subsequently, there have been several 
theoretical studies on odd- frequency pairings [U fl IS H • 
Fuseya, et al. have discussed the possibility of "spin- 
singlet p-wave" pairing near the antiferromagnetic (AF) 
quantum critical point and in the coexisting region with 
AF by using phenomenological pairing interactions in the 
context of the superconducting state in Ce compounds 
0, HI- However, a microscopic mechanism to realize odd- 
frequency pairing has not been established yet, since it is 
not easy to imagine superconductors without equal time 
pairing. 

Although realization of the superconducting state, i.e., 
pair potential in the bulk system, has not been clari- 
fied, odd-frequency pairing correlation, i.e., pair ampli- 
tude, is revealed to be generated in inhomogeneous sys- 
tems where the bulk superconductor is a conventional 
even- frequency one. In ferromagnet/superconductor het- 
erostructures with inhomogeneous magnetization, the 
odd-frequency pairing state was first proposed in Ref. 
|9| and various aspects of this state were subsequently in- 
tensively studied [l(J. Furthermore, the ubiquitous pres- 
ence of odd-frequency pairing amplitude in non-uniform 



superconducting systems has been revealed even if the 
bulk superconductor is a conventional even-frequency one 
[ 1 ll 1 1 2i] . Nowadays, the physics of odd-frequency pairing 
has become an exciting topic for study. 

On the basis of this background, it is an intriguing is- 
sue to clarify whether the odd-frequency bulk supercon- 
ducting state can be derived from a microscopic model 
in a strongly correlated system. It is well known that 
the spin-singlet d 2 .2_ a 2-wave state is favored due to AF 
fluctuations in the Hubbard Model in a square lattice 
fl3l ]. On the other hand, the spin-singlet <i-wave state 
is not favored in a triangular lattice due to geometrical 
frustration [T^|. wherein the nesting condition becomes 
worse. Vojta, et al. discussed the possibility of odd- 
frequency pairing in a triangular lattice @. The merit 
of the odd-frequency pairing is that it can avoid di- 
rect Coulomb repulsion since pairs are formed by two 
electrons at different times by the retardation effect of 
interaction between two electrons. It is highly desir- 
able to clarify favorable conditions which generate odd- 
frequency pairing on the basis of a microscopic calcula- 
tion. 

The purpose of the present study is to explore the odd- 
frequency pairing state in the extended Hubbard model 
on a 2D triangular lattice within the random phase ap- 
proximation (RPA), wherein off-site Coulomb interac- 
tion V as well as on-site Coulomb interaction U exist. 
Near the spin-density- wave (SDW) states, odd- frequency 
pairing becomes dominant. In the presence of V, even- 
frequency spin-singlet d x 2_ y 2-weLve pairing, in which two 
electrons on neighboring sites form pairs, is unfavorable 
since V suppresses the attractive interaction between 
them. By changing the lattice structure from triangular 
to square, even-frequency spin-singlet d x 2_ y 2-w&ve pair- 
ing becomes dominant. Based on these results, preferable 
conditions for the realization of the odd-frequency pair- 
ing are i) geometrical frustration, ii) the presence of V, 
and iii) sufficiently large amount of U which induces the 
SDW order. 
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We consider the extended Hubbard model at half- 
filling on a two-dimensional triangular lattice, where on- 
site repulsive interaction U and off-site interaction be- 
tween nearest neighbor sites V are included. The Hamil- 
tonian is given as 

H = £ (*ij c L c J> + h - c 

* (id) 

where Ci a (') is an annihilation (a creation) operator of 
electron, m a = cl a Ci a and n, = + n^. Within the 
RPA, the spin susceptibility \s and the charge suscepti- 
bility Xc are given as 



Xs(q) 

Xc(q) 



xo(q) 



xo(q) 



l + (U + 2V(q)) X o(q)' 



(2) 
(3) 



where q = (q,ioj m ) = (q,2rmrTi). Here xa is the irre- 
ducible susceptibility given by 



xo(q) 



T 

'N 



(4) 



where k = (k, ie n ) = (k, (2n — l)nTi) and G(k) = (ie n — 
£k + /•*) ■ In the above, n and m are integers, T denotes 
temperature and N denotes the number of sites. In a 
triangular lattice, the band dispersion and the off-site 
Coulomb interaction V(q) are given as 

£fc = —2ti [cos(k x ) + cos(fcz)] — 2t 2 cos(k x + k y ), (5) 
V(q) = 2V{cos(q x ) + cos(q y ) + cos(q x + q y )). (6) 

The effective pairing interactions for the spin-singlet and 
spin-triplet channels due to spin and charge fluctuations 
are given by 

vi ff {q) = u + v(q) + lu 2 Xs (q) 



-(U + 2V{q)?Xc{q), (7) 



V: ff (q) = V(q) - -U 2 Xs (q) 



-(U + 2V{q)) 2 x c {q). 



(8) 



In order to study the pairing instabilities, we solve the 
linearized Eliashberg equation with the power method, 

AA(fc) = -^y e// (fc-fc')G(fc')G(-fc')A(A : '). (9) 



The superconducting transition temperature T c is deter- 
mined by the condition A = 1. We discuss the pairing in- 
stabilities by comparing the magnitude of A for each pair- 
ing symmetry. For the pairing function A(fc), we consider 



the following four symmetries: (i) even-frequency spin- 
singlet even-parity (ESE), (ii) even- frequency spin-triplet 
odd-parity (ETO), (iii) odd-frequency spin-singlet odd- 
parity (OSO), (iv) odd-frequency spin-triplet even-parity 
(OTE). In the present calculation, we choose t\ = 1 for 
the unit of energy. We set the number of electrons as 
n = 1. We take 72 x 72 fc-meshes and 1024 Matsubara 
frequencies. 

First, we focus on the isotropic triangular lattice with 
t\ = *2- In this case, the six irreducible representations 
of group is possible for the pairing function. In this 
paper, we use the notation s- d-, i-, p-, f\- and /2-wave 
for Ai, E 2 , A 2 , Ei, Bi and B 2 , respectively. In Fig. 
[TJ we show the T-dependence of A at half-filling with 
U/ti = 4.0 and V/ti = 0. For each pairing symmetry 
(ESE, ETO, OSO and OTE), we show the result with 
the largest eigenvalue A. For the OSO pairing, the ob- 
tained values of A for the p-wave and /i-wave are very 
close. Thus, we only show the result for the /i-wave. 
It is remarkable that the value of A for OTE pairing 
is the largest for a wide range of temperatures. This 
value has non-monotonic temperature dependence, where 
it reaches unity and then decreases with a decrease in T 
for T/h < 0.1. 



OTE (i-wave) 
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FIG. 1: T-dependence of A at U/h = 4.0, V/ti 
regular triangular lattice. 
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In order to understand the generation of the OTE 
state, we analyze the momentum and frequency depen- 
dence of the pairing function. For this purpose, we trans- 
form Eq. © 



T 



A = -^E/( fc ^')]/[Ei A wi 2 ]- 



(10) 



fe,fe' 



with f(k,k') = V eff (k - k')G(k')G{-k')A(k)A(k'). 
First, we focus on the ESE pairing. Note that there is no 
sign change in A(fc) as a function of the Matsubara fre- 
quency for even-frequency pairing. For spin-singlet pair- 
ing, the pairing interaction V e t / (<?) > is satisfied when 
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the spin fluctuations are dominant. To obtain a large pos- 
itive A, the pairing function A(fc) should change sign by 
the shift of Q ~ (0.47T, OAir), where V e ff(Q,iu> m ) takes 
the maximum value. Therefore, nodes are formed in the 
pairing function consistent with d 2 .2_ 2/ 2-wave pairing. 

Next, we consider the OTE (s-wave). For spin-triplet 
pairing, V e ff(q) always becomes negative. Therefore, the 
pairing function without nodes (s-wave) becomes domi- 
nant. In Fig. [2l we show the Matsubara frequency de- 
pendence of A(k ,ie n ) and V e ff(Q,iu n ) for the OTE 
(s-wave) with U/t 1 = 3.5, T/ti = 0.01. Here, k = 
(7r/2,7r/2) denotes the wave vectors, where A(ko,ie n ) 
has the maximum value. To understand the frequency 
dependence, we ignore the momentum dependence of 
V e ff(q) and A (A;). We introduce the half- value width 
of V e ff(u) m ), r, and the position of peak 7 in A(e„). In 
the summation of the numerator on the right hand side 
of Eq. CG|, V e f f (e n - £„/) for | 

£n ~ £«' I < r mainly 
contributes to A. Since A(e n ) is an odd-frequency of 
E n , f(k,k') with e n e n > > enhances A, while that with 
e n e n i < suppresses A. In order to increase the value of 
A, 7 should be larger than T because the negative parts of 
f(k,k') are reduced with the decrease of \V e ff(ui m )\ for 
cu m > T. On the other hand, G(k')G{-k') = 1/(4 +e|,) 
has a sufficient value only for \e n ' \ ;$ T/t±. If the mag- 
nitude of 7 is larger than T/t\, the positive contribution 
to enhance the magnitude of A becomes small due to the 
small amplitude of A(k')G(k')G(-k'). Therefore, the 
value of A for the OTE pairing decreases at low temper- 
atures with T/ti < 0.1. As a result, a large magnitude 
of A is realized for T jt\ > 7 > T. 
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FIG. 2: Matsubara- frequency dependence of A(ko,ie n ) and 
V e ff(Q,iu n ) for odd- frequency spin-triplet s-wave pairing for 
U/ti = 3.5 and T/ti = 0.01 . 



In Fig. [31 we show the F-dependence of A for the OTE 
(s-wave) and ESE (d-wave) at U/h = 4.0 and T/h = 0.1. 
The value of A for the ESE decreases as V increases. This 
is due to the development of charge fluctuations Xc- For 
the spin-singlet pairing, Xc competes with \s since the 
sign of the coefficients of Xs is opposite to that of Xc in the 



effective interactions. The first order of V, V(q) in Eq. 
(|7|), also competes with Xs- On the other hand, the value 
of A for the OTE slightly increases since the effective 
interaction for the spin-triplet pairing is enhanced by V. 
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FIG. 3: V-dependence of A at U/ti = 4.0 and T/ti = 0.1. 



In the above, we have shown that the odd-frequency 
pairing is favored near the SDW state on a triangular lat- 
tice. In the following, we clarify how the above results are 
influenced by the lattice structure. In Fig. [H we show the 
i2-dependence of A for OTE (s-wave) and ESE (d x 2_ y 2- 
wave), where t% is a transfer integral on the diagonal. 
This lattice corresponds to the square lattice in the case 
of t2/t\ = 0, and it corresponds to the triangular lattice 
in the case of t2/*i = 1- The critical point of the SDW 
state is different in each lattice, since the Xoil) value 
obtained is different. Here, we introduce the Stoner fac- 
tor, sf — Uxmax, which gives the closeness to the SDW 
state, where Xmax is the maximum of Xo(<?)- s f — 1 cor- 
responds to the SDW critical point, and sf is less than 
unity in the paramagnetic phase. We fixed the Stoner 
factor to be sf = 0.92 for varying the values of tijt\ for 
the purpose of keeping the distance from the SDW crit- 
ical point. The values of A for both OTE (s-wave) and 
ESE (c? a .2_ y 2-wave) increase with ti for £2/^1 < 0.7, since 
U, which is determined by U = sf /xmax, increases with 
t 2 - The value of A for ESE (d x 2_ y 2-w&ve) decreases with 
t 2 in the region t%Jt\ > 0.8, while the magnitude of U 
used in the actual calculation increases with t 2 . This is 
because geometrical frustration becomes significant with 
increasing t% and the resulting nesting condition becomes 
worse. On the other hand, A for OTE (s-wave) contin- 
ues to increase with t 2 in the region t%jt\ > 0.8, since 
the frustration does not affect the instability of s-wave 
pairings. 

Finally, we show that odd-frequency pairing can be 
realized in a frustrated square lattice with next nearest 
neighbor hopping (t 2 ), where the band dispersion is given 
as = — 2ii(cos(fc E ) + cos(fc x )) — 4i 2 cos(k x ) cos(k y ). In 
Fig. O the T-dependence of A is shown for t%jt\ = —0.6, 
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FIG. 4: £ 2 -dependence of A at sf = 0.92 and T/ti = 0.2. 



U/ti = 2.88. It can be shown that OTE (s-wave) is dom- 
inant at high temperatures as in the case of the triangular 
lattice. The AF fluctuation is suppressed due to the ge- 
ometrical frustration induced by large £2 ■ Therefore, the 
value of A for ESE d x 2_ y 2--wave pairing decreases and be- 
comes smaller than that of the OTE pairing. Note that 
V enhances the value of A of the d xy -~w&ve pairings in 
the square lattice while it suppresses that of the d x 2_ y 2- 
wave pairing. However, the value of A for the d xy -wave 
pairings remains smaller than that of the d x 2_ y 2-wave. 
Thus, we suggest that odd-frequency pairing appears in 
geometrical frustrated systems. 
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In conclusion, we have investigated whether odd- 
frequency pairing states can be realized in the extended 
Hubbard model at half-filling on a 2D triangular lattice 
using the RPA. The odd-frequency spin-triplet s-wave 
pairing is dominant for a wide temperature regime. The 
instability of the spin-triplet s-wave state is enhanced by 
the off-site Coulomb interaction V because the effective 
interaction for the triplet channel is enhanced by V . On 
the other hand, since the effective interaction for the spin- 
singlet channel is reduced by V, the even-frequency spin- 
singlet d-wave state is suppressed. It is an interesting 
future problem to study the realization of odd-frequency 
pairs in the presence of a magnetic field or modulation 
of the pair potential. In addition, clarifying the super- 
conducting properties of odd-frequency superconductors 
is also an interesting issue. Since phase sensitive probes, 
e.g., tunneling and Josephson effects, are crucial to iden- 
tify the pairing symmetry [l], [l5| in unconventional su- 
perconductors, similar studies on odd-frequency super- 
conductors will become important (l6l . Il7j . 
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